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Abstract 

Chaotic walking of cold atoms in a tilted optical lattice, created by two counter propagating running 
waves with an additional external field, is demonstrated theoretically and numerically in the semi- 
classical and Hamiltonian approximations. The effect consists in random-like changing the direction 
of atomic motion in a rigid lattice under the influence of a constant force due to a specific behavior 
of the atomic dipole-moment component that changes abruptly in a random-like manner while atoms 
cross standing-wave nodes. Chaotic walking generates a fractal-like scattering of atoms that manifests 
itself in a self-similar structure of the scattering function in the atom-field detuning, position and 
momentum spaces. The probability distribution function of the scattering time is shown to decay in 
a non-exponential way with a power-law tail. 
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1 Introduction 



The mechanical action of light upon neutral atoms placed in a laser standing wave is at the heart 
of laser cooling, trapping, and Bose-Einstein condensation [lj. Numerous applications of the mechan- 
ical action of light include isotope separation, atomic interferometry, atomic lithography and epitaxy, 
atomic-beam deflection and splitting, manipulating translational and internal atomic states, measure- 
ment of atomic positions, etc. Atoms and ions in an optical lattice, formed by a laser standing wave, 
are perspective objects for implementation of quantum information processing and quantum computing. 
Advances in cooling and trapping of atoms, tailoring optical potentials of a desired form and dimension, 
controlling the level of dissipation and noise are now enabling the direct experiments with single atoms 
to study fundamental principles of quantum physics, quantum chaos, decoherence, and quantum-classical 
correspondence. 

Nonlinear dynamics of cold atoms in optical lattices is a fastly growing branch of atomic physics. 
There are a number of theoretical works and impressive experiments on quantum chaos, dynamical 
localization, chaos-assisted tunneling, Levy flights, etc. (for reviews see [HE]). To suppress spontaneous 
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emission and provide a coherent quantum dynamics one usually works far from the optical resonance. 
Adiabatic elimination of the excited state amplitude leads to an effective Hamiltonian for the center- 
of-mass motion, whose 3/2 degree-of- freedom classical analogue has a mixed phase space with regular 
islands embedded in a chaotic sea. New possibilities are opened if one works near the optical resonance 
and take the internal atomic dynamics into account. A single atom in a standing- wave laser field may be 
semiclassically treated as a nonlinear dynamical system with coupled internal (electronic) and external 
(mechanical) degrees of freedom [3H6]. I n the semiclassical and Hamiltonian limits (when one treats 
atoms as point-like particles and neglects spontaneous emission and other losses of energy), a number 
of nonlinear dynamical effects have been analytically and numerically demonstrated with this system: 
chaotic Rabi oscillations [IHS], Hamiltonian chaotic atomic transport and dynamical fractals [THlOj. Levy 
flights and anomalous diffusion [6lllltfT2]. These effects are caused by local instability of the CM motion 
in a laser field. A set of atomic trajectories under certain conditions becomes exponentially sensitive to 
small variations in initial quantum internal and classical external states or /and in the control parameters, 
mainly, the atom-laser detuning. Hamiltonian evolution is a smooth process that is well described in 
a semiclassical approximation by the coupled Hamilton-Schrodinger equations. A detailed theory of 
Hamiltonian and dissipative chaotic transport of atoms in a laser standing wave has been developed in 
Refs. P3J] and [HUH], respectively. 

Additional possibilities to manipulate the atomic transport are created by applying an external force 
to the standing-wave optical potential. It is obvious that for cold atoms in a vertical optical lattice it is 
necessary to account for the Earth's acceleration. It is possible as well to create horizontal accelerated 
optical lattices by adding a constant force whose magnitude along the optical axis can be easily varied. 
The problem of atomic motion in a tilted optical potential is closely related to the old problem of electron 
motion in a a periodic crystal with dc or ac forces applied. The analogue of well known Bloch oscillations 
with cold atoms has been experimentally found in Ref. [15|J16|. 

In the present paper we apply the ideas and methods, elaborated in the field of nonlinear dynamics of 
cold atoms, to study theoretically and numerically motion of point-like atoms in a tilted optical lattice. 
It will be shown that varying only one parameter, the detuning between the frequencies of a working 
atomic transition and the laser field, one can explore a variety of regimes of atom motion, including 
chaotic walking, dynamical fractals and chaotic scattering. 



2 Chaotic and regular regimes of motion of atoms in a tilted potential 



In the one-dimensional case, the Hamiltonian of a two-level atom in a standing-wave laser field and 
an additional external field can be written in the frame rotating with the laser frequency ujf as follows: 

P 2 h 

H = - + -{oJ a - 0J f )a z - nn (<r_ + &+) cosk f X + FX, (1) 

2m a 2 

where fr± x are the Pauli operators for the internal atomic degrees of freedom, X and P are the classical 
atomic position and momentum, uj a and f2o are the atomic transition and maximal Rabi frequencies, 
respectively. F stands for the static force induced by external field. 

In the semiclassical approximation, where the transversal atomic momentum p is supposed to be, in 
average, much larger than the photon one Hkf, atom with quantized internal dynamics is treated as a 



2 



point-like particle to be described by the Bloch-Hamilton equations of motion without relaxation terms 

x = uj r p, p = — usinx — K, u = Av, 

(2) 

v = —Au + 2z cos x, z = —2v cos x, 

where u and v are synchronized (with the laser field) and quadrature components of the atomic electric 
dipole moment, respectively, and z is the atomic population inversion. Equations (|2|) are written in the 
dimensionless form with x = kfX and p = P/hkf to be classical atomic center-of-mass position and 
momentum, respectively. Dot denotes differentiation with respect to the dimensionless time r = Qt. The 
set ([2]) has the three control parameters 

uj r = hk 2 /m a £l, A = (uif — uj a )/n, n = F/hkftt, (3) 

which are the normalized recoil frequency, co r , atom-field detuning, A, and applied force k. The system 
has two integrals of motion, namely the total energy 

ti = —p +KX — UCOSX — —z, (4J 

and the length of the Bloch vector u 2 + v 2 + z 2 = 1. 

The external force is directed in the negative direction of the optical axis x. So, if the initial atomic 
momentum, po, is chosen to be in the negative direction, the force will simply accelerate the corresponding 
atoms. If po > 0, then one may expect much more complicated motion. 

Equations ([2]) constitute a nonlinear Hamiltonian autonomous system with two and half degrees of 
freedom. Owing to two integrals of motion, phase points move on a three-dimensional hypersurface 
with a given energy value H. In general, motion in a three-dimensional phase space in characterized 
by a positive Lyapunov exponent, a negative exponent equal in magnitude to the positive one, and zero 
exponent. The maximal Lyapunov exponent characterizes the mean rate of the exponential divergence 
of initially close trajectories and serves as a quantitative measure of dynamical chaos in the system. 
Because of a transient character of chaos, we have computed the finite-time Lyapunov exponent A by the 
method developed in Refs. [17^118] . It has been found that at the fixed value, uj r = 10 -3 , of the recoil 
frequency A is positive in the following ranges of values of the other control parameters: the detuning 
—0.5 < A < 0.5 and the force —0.25 < k < 0.6. Therefore, we expect chaotic motion of atoms with the 
parameter's values in those ranges. 

In numerical experiments throughout the paper we suppose that two-level atoms are initially prepared 
in the ground states, uq = v o = 0, zq = —1, with xq = and fixed values of the two control parameters, 
the normalized recoil frequency, uv = 10 -3 , and the external force n = 0.01. The atom-field detuning, A, 
can be changed in a wide range. It will be shown in this paper that atoms may perform chaotic walking, 
the new type of motion in absolutely deterministic environment where atoms can change the direction 
of motion alternating between flying through the standing wave and being trapped in its potential wells. 
We would like to stress that the local instability produces chaotic center-of-mass motion in a rigid optical 
lattice without any modulation of its parameters. 

In Fig. [1] we illustrate different regimes of the center-of-mass atomic motion along the optical axis with 
the initial atomic momentum chosen to be po = 10. It is simply the motion on the phase plane x — p. A 
typical picture with chaotic walking is shown in Fig. [Ih. with the value of the detuning, A = 0.15, at which 
the maximal Lyapunov exponent, A, is positive. The atom starts to move in the positive x-direction, 



3 



changes soon the direction of motion a few times, acquiring irregularly positive and negative values of 
the momentum, and suddenly begins to move in the positive x-direction for a comparatively long time. 
Then it is decelerated, turns back and flies in the negative direction. After that it changes the direction 
of motion many times demonstrating what we call "chaotic walking" . 

For comparison, we show in Fig. [lb the phase-plane motion with the larger detuning A = 1 (and 
at the same other conditions) at which the maximal Lyapunov exponent is not positive. The atom 
moves initially in the positive x-direction, accelerating and decelerating alternatively. Soon it changes 
the direction of motion and moves permanently in the negative direction. The motion is regular with 
a slight modulation of the momentum. The motion at exact resonance, A = 0, is even more simple 
(Fig. Eh). 

What is the ultimate reason of chaotic walking? For an optical lattice without an external force, it has 
been found in Ref. [10] that instability is caused by the specific behavior of the Bloch-vector component of 
a moving atom, u, whose shallow oscillations between the standing-wave nodes are interrupted by sudden 
jumps with different amplitudes while atom crosses each node of the wave. It looks like a random like 
shots happened in a fully deterministic environment. The reason of chaotic walking in a tilted potential 
is the same. It follows from the second equation in the set ([2]) that those jumps result in sudden changes 
of the atomic momentum while crossing nodes. If the value of the atomic energy is close to the separatrix 
one, the atom after the corresponding jump-like change in p can either overcome the potential barrier 
and leave a potential well or it will be trapped by the well, or it will move as before crossing the node. 
The evolution of all the Bloch components in the regime of chaotic walking is shown in Fig. [2j For 
comparison, we show in Figs. [3] and [4] their evolution in the regular regimes, far off the resonance and 
at exact resonance, respectively. 

3 Fractal scattering of atoms in a tilted potential 

Different types of fractal- like structures may arise in chaotic Hamiltonian systems (see reviews [19|I20| ) . 
It is known from many studies in celestial mechanics [21], fluid dynamics [221123], atomic physics [711101 
[]2J[25], cavity quantum electrodynamics [8,9], underwater acoustics [21] and other disciplines [26] that 
under certain conditions the motion inside an interaction region may have features that are typical for 
dynamical chaos, (homoclinic and heteroclinic tangles, fractals, strange invariant sets, positive finite-time 
Lyapunov exponents, etc.) although the particle's trajectories are not chaotic in a rigorous sense because 
chaos is defined as an irregular motion over infinite time. 

Let us place atoms one by one at the point xq = with the same value of the initial momentum 
Po = 10 but change slightly the value of the detuning A. All the other initial conditions and the control 
parameters are supposed to be the same for all the atoms. We fix the time moment T when each 
atom crosses the point x = moving in the negative direction. The exit time function T(A) in Fig. [5] 
demonstrates the complicated structure with smooth intervals alternating with wildly oscillating peaks 
that cannot be resolved in principle, no matter how large the magnification factor. The panels (b) and 
(c) in Fig. [5] are successive 50 times magnifications of the detuning intervals shown in the panel (a). 
Further magnifications reveal a self-similar fractal-like structure that is typical for Hamiltonian systems 
with chaotic scattering. The exit time T increases with increasing the magnification factor. The same 
picture is observed when computing the exit time function in the position and momentum spaces. It is 
a clear demonstration of a fractal-like behavior of chaotically walking atoms. 
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It is established in theory of one and half degree-of-freedom systems that transient Hamiltonian 
chaos in the interaction region occurs due to existence of, at least, one non-attractive chaotic invariant 
set consisting of an infinite number of localized unstable periodic orbits and aperiodic orbits. This set 
possesses stable and unstable manifolds extending into the regions of regular motion. The particles with 
the initial positions close to the stable manifold follow the chaotic-set trajectories for a comparatively 
long time, then deviate from them, and leave the interaction region along the unstable manifold. In a 
typical Hamiltonian system there exists an infinite number of trajectories of zero measure with infinite 
exit time which belong to that chaotic invariant set. Our system with two and half degrees of freedom is a 
much more complicated one, and it is practically impossible to reveal the corresponding chaotic invariant 
set with its stable and unstable manifolds. However, the mechanism of chaotic scattering and fractal-like 
structures should be the same. 

The statistics of exit times T is shown in Fig. [6] in a semilogarithmic and logarithmic scales. The 
probability distribution function (PDF) in this figure gives the probability for an atom to have a given 
value of T. The bold straight line in Fig. [6^ implies that the PDF is exponential in its middle part, 
P ~ exp(— aT), with the exponent a = —0.000270722. However, the tail of the PDF is not exponential. 
To prove that we plot the function in the logarithmic scale in Fig. [BJa and compute the slope at the tail. 
It has no sense to calculate the slope at the very tail because of a small number of events with very 
large values of T. The bold straight line implies that the PDF is a power-law one, P ~ T~ 7 , with the 
coefficient 7 = —2.53086. It is interesting that the slope at the PDF tail around the value —2.5 is rather 
typical for many chaotic Hamiltonian systems [27, 28]. The reason of that is unclear. 

In hyperbolic chaotic systems the PDFs should decay exponentially because the phase space of such 
systems is homogeneous, and all the trajectories are unstable. It is not the case even with one and 
half degree-of-freedom systems with inhomogeneous phase space, where exist so-called stability islands 
embedded in a stochastic sea, because the borders of those islands are "sticky" . It means that a typical 
chaotic trajectory, wandering in the stochastic sea, approaches the island's borders and "stick" to them 
for a long time. By that reason, the corresponding PDFs are not exponential but power-law ones at their 
tails. PDFs with power-law decay imply that the corresponding quantity, the exit time in our case, is 
scale invariant i.e., there is no a single dominant scale in the process. Geometrically it means that chaotic 
trajectories for such a process are self-similar. 

4 Conclusion 

It is shown that point-like atoms in a tilted optical potential with a constant external force applied 
can move chaotically changing the direction of motion in a random-like way. The existence of chaos 
is confirmed by direct computation of the maximal finite-time Lyapunov exponent of the equations of 
motion that is shown to be positive in a range of the atom-laser detuning and the applied-force strength. 
The ultimate reason of chaotic walking is the specific behavior of the Bloch- vector component of a moving 
atom, u, whose shallow oscillations between the standing-wave nodes are interrupted by sudden jumps 
with different amplitudes while atom crosses each node of the wave. It is demonstrated numerically that 
such a behavior arises exactly at those values of the detuning for which the Lyapunov exponent is positive 
and atoms move chaotically. We illustrate different regimes of the center-of-mass motion simply varying 
the detuning. It is an easy way to manipulate the atomic transport in tilted optical lattices. 

Treating motion of atoms in a tilted optical lattice as a scattering problem, we show that the scattering 
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Figure 1: Motion of a cold atom in a deterministic tilted optical lattice as it looks on the phase plane 
x — p. (a) Chaotic walking at A = 0.15, k = 0.01, u r = 10~ 3 . (b) Regular motion at A = 1 and with the 
same other conditions, (c) Regular motion at the resonance, A = 0, with the same other conditions. 




Figure 2: Evolution of the atomic Bloch components in the regime of chaotic walking (A = 0.15). 

of atoms under conditions of chaotic walking is chaotic and typical for Hamiltonian systems. Fixing 
the time moment T when atoms with slightly different values of the detuning, momentum or initial 
position cross a fixed point (x = 0), we show that the corresponding scattering functions demonstrate 
the complicated structure that cannot be resolved in principle, no matter how large the magnification 
factor. Owing to that the probability to have a given value of T is not exponential but decays at its tail 
by a power law. 
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Figure 5: Atomic dynamical fractal. Self-similar dependence of the exit time, T, with given initial 
position, xq = and momentum po = 10, on the detuning. The successive magnifications are shown. 
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